Onset of Irreversibility in Cyclic Shear of Granular Packings 
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We investigate the onset of irreversibility in a dense granular medium subjected to cyclic shear 
in a split-bottom geometry. To probe the micro and mesoscale, we image bead trajectories in 3D 
throughout a series of shear strain oscillations. Though beads lose and regain contact with neighbors 
during a cycle, the global topology of the contact network exhibits reversible properties for small 
oscillation amplitudes. With increasing reversal amplitude a transition to an irreversible diffusive 
regime occurs. 

PACS numbers: 45.70.-n 47.57.Gc 64.60.aq 64.60.ah 

The transition from reversible to irreversible dynam- 
ics in systems of many interacting particles is of fun- 
damental importance to many-body physics. Granular 
materials under repeated shear deformations constitute 
excellent experimental systems with which to study such 
reversibility for out-of-equilibrium systems. Since grains 
are macroscopic in size, thermal fluctuations are negligi- 
ble [H, with particle trajectories determined by contact 
forces between neighboring grains. In ground breaking 
work on low concentration granular suspensions, Corte 
and coworkers observed that particle collisions tend to 
disrupt the inherent reversibility of low Reynolds num- 
ber fluid flow, and that reversibility is lost at strains and 
densities that allow a particle to collide with two or more 
neighbors Q • As a consequence, the onset of irreversibil- 
ity decreases approximately as the square of particle den- 
sity i. 

For dense particle systems close to jamming, the 
threshold for exact spatial reversibility would thus be- 
come extremely small. However, dense configurations of 
particles have additional constraints on their positions: 
the neighborhood of nearby grains in the local environ- 
ment of the particle confine its position, providing a new 
reference for the motion of the particle. This perspec- 
tive opens up the possibility of an alternative form of 
reversibility with respect to the local neighborhood of 
particles. Under shear strain, neighbors are lost, and 
correlations are important: rearrangements in dense con- 
figurations are highly collective and are known to exhibit 
a range of dynamical features, from local reversibility 
as found in foams \^ to larger scale dynamical hetero- 
geneities as found in 2D bidisperse cylinders 

In this article, we describe experiments that explore 
strains for which these rearrangements first emerge at 
scales larger than a grain diameter. Since we probe ex- 
tremely slow strain rates, the granular temperature is 
virtually zero and rearrangements are due to plastic de- 
formations as opposed to random motion. We character- 
ize the local environment around particles, and the plas- 
tic deformations, through the analysis of the network of 



nearby particles. In our experiments, the full three di- 
mensional motion of beads within a split-bottom shear 
cell undergo cyclic shear. Dense granular systems under 
cyclic shear in two and three dimensions have been stud- 
ied previously 0, [§"111: but the trajectories of all beads 
in the bulk have not been probed, and collective bead 
dynamics and reversibility are not well understood. 

We find that Mean Square Displacements (MSD) of 
beads indicate a crossover between two different behav- 
iors and that a characteristic strain amplitude is associ- 
ated with this crossover. In order to relate this change in 
behavior to changes in the neighbor network, we analyze 
neighbor changing events through the use of the Broken 
Link Network [12l | , which assigns a link to pairs of parti- 
cles that then later move apart, under shear strain. For 
a strain amplitudes below the aforementioned character- 
istic amplitude, the largest interconnected component of 
the BL network nearly disappears after shear reversal, 
while above the characteristic amplitude a giant compo- 
nent appears to emerge and survive the shear reversal. 
We posit a novel type of reversibility/irreversibility tran- 
sition under cyclic shear strain: one that distinguishes 
between a regime where rearrangements are essentially 
local and random, never leading to a qualitative change 
in the topology of the BL network, and a regime where 
the topology is modified significantly by the formation 
of a giant component. We refer to this phenomenon as 
'topological irreversibility,' by virtue of the fact that the 
topology of the BL network is being irreversibly altered 
by the formation of a giant component. 

Setup — To study motion of all beads, we built a 



modified version of the split-bottom geometry [13|, |14| 
that allows for 3D imaging. In this geometry, shear is 
applied from below by rotating a disk mounted flush with 
the bottom. A shear zone is then emanating from the 
edge of a rotating disk — much of the shearzone is far 
from boundaries, the shear zone is wide and smooth, and 
nothing obstructs its view, making it ideally suited for 
3D observations of grain flows with smooth gradients. 
We use a square box with a 15 cm by 15 cm bottom 
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the fraction of grains undergoing differential flow. 

The angular velocity flow field is quite heterogeneous in 
the case of constant shear. The differential flow, as shown 
in [1^ is greatest in a region above the shearing disk. We 



FIG. 1: (a) Split-bottom shear cell with acrylic beads in index 
matched, dyed fluid. A laser sheet illuminates cross-sections 
for imaging, (b) As shown in this reconstruction, the po- 
sition of all beads is found by scanning the laser sheet sheet 
across the whole container and automatic particle recognition. 

(c) -(d) A typical trajectory of a bead under oscillatory shear 
strain with a reversal amplitude of 40° from (c) above and 

(d) the side (cube sides are 2 bead diameters, color changes 
from blue to red color with increasing time). 



and a 9 cm diameter disk at the bottom. In order to 
provide traction for the beads above, we attach a sparse 
coating of 3 mm beads to the bottom of the container 
and disk surface. We fill the container with 3/16" diam- 
eter acrylic beads. We immerse the beads in an index 
matched fluid (Triton X-100). A fluorescent dye (Nile 
Blue 690 Perchlorate) is added to the fluid, along with 
a small amount of HCl (1.0 mL per 1.0 L of fluid) to 
stabilize the mixture. 

We use the 3D imaging technique described in [l5j and 
illustrated in Fig. [1] but with an improved bead finding 
algorithm. We use a 3D version of the convolution de- 
scribed in [3l to find the position of each bead to better 
than half a pixel, or 4% of a bead diameter. We then 
track the motion of individual beads with a tracking al- 
gorithm developed by Grier et al. 17[ to find trajectories 
such as the one in Fig. [TJ:-d. With the full trajectories 
for about 98% of beads in our system outside the bottom 
layer, we can now analyze how beads move with respect 
to their neighbors. 

We choose to fill the container to a filling height H 
equal to the radius of the shearing disk Rg . This regime, 
which has only been approached experimentally in a lim- 
ited number of studies 12, 18, 31, is between two pre- 
dicted regimes wherein the shear bands in a constant 
strain angular velocity profile transition from annular re- 
gions that widen above the disk edge {H/Rg < 0.7) to 
a dome-shaped region above the disk (H/Rg » 1) [20|. 
We choose this intermediate regime in order to maximize 



choose to limit our analysis to this region of maximum 
relative strain. We find this region by calculating the 
mean angular velocity in annular regions about the axis 
of disk rotation. We disregard axial and radial velocities, 
since the bulk remains stationary on average. We use 
these values to find the approximate differential strain in 
the axial and radial directions. The regions with relative 
strains in the lower 60th percentile are excluded. 

Imaging cannot distinguish whether two hard spheres 
are in contact or just close together; instead we choose 
a maximum separation criterion to identify the nearest 
neighbors of each bead based on a pref erence for tan- 
gential motion of neighboring beads |l2l . [ilj . Based on 
this criterion, we consider beads whose centers are within 
w 1.08 bead diameters of one another to be nearest neigh- 
bors, since beads that are separated by greater distances 
show no such preference for tangential relative motion. 
Thus, the beads have an effective radius R ~ 0.54 bead 
diameters. 

Protocol — To create reproducible initial conditions 
representative of steady shear, we rotate the disk two 
full revolutions at 1 mrad/sec. At this rate of strain, the 
granular flow in our fluid immersed system is indepen- 
dent of strain rate, with the same velocity profile as a 
dry granular flow 2M- then reverse the direction of 
rotation of the disk and move the disk for a strain of Or 
(= 2°, 4°, 10°, 20°, 40°). We reverse the direction again 
and rotate the disk by —0r to its original position, thus 
straining and reversing the system in a cyclic fashion. We 
perform 20 successive cycles in each experiment, taking 
a 3D image of the system after every 2° of strain. 

Mean Square Displacements — The first measure of 
reversibility we employ is to test how far a bead has 
moved after a full cycle. The mean square displacements 
(MSD) of the beads in the shear zone at the end of each 
of the 20 cycles are shown as function of cycle number, 
TT-c, in Fig. [5^. For all but the largest reversal ampli- 
tude 6r, the slope of the MSD decreases with increasing 
cycle number to slopes below 1, signaling sub-diffusive 
behavior. 

To explore whether the system becomes more re- 
versible after several shear cycles, we choose the begin- 
ning frame of a later starting cycle reference 
frame. Fig.[2}3 illustrates that for rig = 10, the MSD is up 
to an order of magnitude smaller, suggesting that bead 
positions in this cyclic shear flow become more reversible 
after several cycles. Strikingly, the MSD for the largest 
reversal amplitude [Or ~ 40°) is considerably larger in 
magnitude. 

To compare the rearrangements after one cycle of am- 
plitude Or with two cycles of amplitude 9^/2 etc, we cal- 
culate the MSD after a fixed absolute strain /S.9 (relative 
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FIG. 2: (a)-(b) Mean Square Displacements (MSDs) at the 
end of each cycle as a function of cycle number for different 
cycle amplitudes, scaled by the effective bead radius R: (a) 
Reference frame at the start of the first cycle, (b) The start of 
the f 0th cycle is chosen as a reference frame, (c) MSD after 
80° of absolute strain vs. reference cycle number. 



to a chosen reference frame) which we fix here at 80°. 
Fig. [5}; shows the MSD at this fixed total strain as func- 
tion of the start cycle Ug. We see that for — 40°, the 
magnitude of the MSD does not depend strongly on the 
choice of reference cycle after a short transient of about 
3 cycles, while the MSD for smaller continue a down- 
ward trend for all cycles investigated. 

In terms of spatial reversibility where all beads return 
to their original coordinates, all three probes of the MSD 
strongly suggest a qualitative change in behavior between 
Qr = 20° and = 40°: for large amplitude, the motion 
is irreversible, and reaches a steady, diffusive state, while 
for small amplitudes, the cyclic shear strain process be- 
comes increasingly reversible (i.e. decreasing MSD mag- 
nitudes) and the system is aging. 

Broken Links — As a second measure of reversibility 
we probe the network of nearest neighbors. We first iden- 
tify all neighbors in a reference frame at the start of a 
cycle. Then we calculate, as function of frame number, 
the fraction of links broken ft, , defined as the fraction of 
neighbor pairs from the reference frame that are no longer 
within the cutoff distance used to define neighbors. The 
approach of 12| has been modified here to include the 
possibility of "healing." A link is healed when a neigh- 
bor pair from the reference frame that has been broken 
later reforms. A healed link does not appear in the bro- 
ken link network. Note that other pairs of beads may 
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FIG. 3: (a) Fraction of links broken, for shear reversal vs. 
2° strain step for the case of dr = 10°. (b) /(, after one cycle 
as a function of the starting cycle n^, for all 9r. (c) ft, as a 
function of AS for representative curves at ns=5, 10 and 15. 



become neighbors in later frames but this is not consid- 
ered in the broken link analysis. 

In Fig. |3K we show the evolution of fh as a function of 
A0, for Or — 10° — other Or show comparable results. 
The top curve shows /& with ns=l. We observe a strong 
decrease in when the strain direction is reversed (at 
10°) — broken links reform after reversal. Additionally, 
fewer additional links are broken in each subsequent cycle 
— the motion becomes less irreversible for later times. 
The ten lower curves in Fig. show when the start 
of a later cycle is chosen as a reference frame, illustrating 
that the recovery of neighbors after a full shear strain 
cycle (at arrow) improves with increasing cycle number. 

In Fig. we show after one cycle as a function 
of the starting cycle n^, which provides a local measure 
of the reversibility of the configuration after one cycle. 
Clearly, the fraction of links broken decreases with in- 
creasing Us (indicating an increase in reversibility) for all 
but the largest reversal amplitude, which appears to level 
off. As a simple local metric, /{, captures what fraction 
of the neighbor network is evolving with strain, a topo- 
logical analog to the "particle activity" of 2]. Indeed, 
the apparent power dependence of the decay on Ug in 
the reversible regime, and a leveling off for the highest 
reversal amplitude, is consistent with the observation of 
a reversible-irreversible transition observed in [5]. How- 
ever, as illustrated in Fig.[3J:, fb cannot identify this local 
reversibility in the neighbor network over the period of 
one cycle - /{, after one cycle increases approximately lin- 
early with cycle amplitude without a clear change in the 
form of the dependence on AO. 

Our characterization of the irreversibility revealed a 
striking difference between the growth of mean square 
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FIG. 4: (a) The peak largest component size versus the 
final largest component size, for each cycle, (a, Inset) Sg 
as a function of A6' over one cycle for Us —5, 10, and 15. 
Peak and end values for Sg are pointed to for the 40° case, 
(b) Sg as a function of the fraction of links broken for 
each reference network used, for 6r of 10° (orange triangles), 
20° (green squares) and 40° (blue circles). 



displacements and fraction of links broken with driving 
amplitude 9r- The MSD barely grows with 6r for small 
driving amplitudes (see in particular Fig. [2t for small 
start cycles, where all data for dr from 4° to 20° strongly 
overlap), but then exhibits a jump for 0r ~ 40°. In con- 
trast, the fraction of links broken /f, grows approximately 
linearly with 9r (see Fig. [St) and so over a single strain 
cycle fb fails to distinguish Or = 40 from the smaller re- 
versal amplitudes. 

We believe that this difference occurs along with the 
percolation of the network of links broken. In this pic- 
ture, a small fraction of links broken leads to small clus- 
ters of failures that upon reversal are mostly healed. Con- 
sequently, after a cycle the beads return to their start- 
ing position, leading to a small MSD. However, once the 
fraction of links broken is sufficiently large, a dominant 
failing cluster forms. Upon reversal, this dominant fail- 



ing cluster is not fully healed, and after one cycle, beads 
end up with different neighbors and at different locations, 
leading to a large MSD. 

To probe this picture, we investigate how links that 
break are connected with each other by observing an as- 
pect of the topology of the broken links network. Un- 
der steady shear, it is known that the size of the largest 
component, i.e, the largest cluster of beads exhibiting 
inter-connected broken links, shows a transition beyond 
a characteristic strain (see [12[). Here we adapt this mea- 
sure to oscillatory shear strain, and follow Sg, which is 
the fraction of beads in the largest component, as func- 
tion of absolute strain A9, with reference frames chosen 
at the beginning of each cycle. 

Fig. |4] provides qualitative support for our picture. In 
the inset of Fig. |4^ we show that for 9r < 20°, the size of 
the largest component Sg initially grows with A9 but then 
shrinks back to almost zero when the motion is reversed 
(links are reformed after reversal), while for 9r = 40°, Sg 
remains substantial after completing a full cycle — hence 
links remain broken. To stress this qualitative difference, 
we collect the data for the last 15 cycles, characterize the 
functional form of Sg{A9) by its peak size at the reversal, 
gpeak^ and the largest component size at the end of the 
cycle, Sg""*, and make a scatter plot of g^^"'"' versus s™'' 
(Fig. 13^). This reveals a clear separation of the 10° and 
20° data which form peak component sizes that can be 
almost fully reduced to zero, from the largest amplitude 
40° where broken links interconnect a significant part of 
the system at the end of a cycle. 

In Fig. |4)d we show scatter plots of Sg versus fb for 
three cycle amplitudes, which show that for each cycle 
amplitude, the data collapses onto a master curve quan- 
titatively similar to that of For the largest driving 
amplitude, a regime emerges where Sg has a dependence 
on fb that is consistent with a power law — a collective, 
percolation effect. 

Discussion — We have studied a granular system un- 
der shear reversal at microscopic and mesoscopic scales. 
Such dense systems are never found to be strictly re- 
versible, but we have shown that irreversibility comes in 
two distinct flavors: spatial and topological. For small 
driving amplitudes, neighbor links are broken but most 
reform, there are no large clusters, and the mean square 
displacements are small, and importantly, the cycle-to- 
cycle mean square displacements decrease with cycle 
number: the system is aging. However, for large driv- 
ing amplitudes, a substantial fraction of neighbor links 
do not reform upon reversal, and the cycle-to-cycle mean 
square displacements become independent of cycle num- 
ber: here the motion becomes diffusive. Our data is con- 
sistent with the hypothesis that this strong irreversibility 
at large cycle amplitudes is connected to the growth of a 
giant component, which signifies a collective breaking of 
contacts involving a number of beads that is proportional 
to the system size. We further showed that the growth 



5 



of the largest component is consistent with a power law 
dependence on fh for the largest reversal amplitude. 

We note that our new topological measures of irre- 
versibility — though less stringent than exact reversal 
symmetry of the equations of motion — correlate well 
with a strong increase in MSD with cycle amplitude 
(Fig. [2];), and thus are likely relevant for understanding 
dynamics of dense granular flows. 

One important future direction of study will be to con- 
nect the onset of irreversibility in bead configuration to 
the emergence of dynamical heterogeneities. It is impor- 
tant to note that this work has probed in detail time 
scales much shorter than those at which dynamical het- 
erogeneities are observed in, for instance, [7|. As a result, 
our observations are likely related to the formation of the 
building blocks of dynamical heterogeneities, whose dy- 
namics may unfold at much larger time scales. 

We thank Joost Weijs for his contributions to our par- 
ticle extraction routine, William Derek Updegraff for his 
assistance in experimentation, as well as Mark Hererra 
for his assistance in analysis. This work was supported 
by NSF grant DMR0907146. KR was supported by the 
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